Tests of mode-coupling theory in two dimensions 
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We analyze the glassy dynamics of a binary mixtures of hard disks in two dimensions. Predictions 
of the Mode-Coupling theory (MCT) are tested with extensive Brownian dynamics simulations. 
Measuring the collective particle density correlation functions in the vicinity of the glass transition 
we verify four predicted mixing effects. For instance, for large size disparities, adding a small amount 
of small particles at fixed packing fraction leads to a speed up in the long time dynamics, while at 
small size disparity it leads to a slowing down. Qualitative features of the non-ergodicity parameters 
and the /3-relaxation which both depend in a non-trivial way on the mixing ratio are found in the 
simulated correlators. Studying one system in detail we are able to determine its ideal MCT glass 
transition point as ip'^ = 0.7948 and test MCT predictions quantitatively. 

PACS numbers: 64.70.P-, 64.70.Q-, 82.70.Dd 



I. INTRODUCTION 

Simulations of simple model systems provide a pow- 
erful and important mean to study the interesting phe- 
nomenon known as glass transition Considering a 
binary mixture of hard disks close to vitrification in a 
computer simulation unveils the opportunity to investi- 
gate a system, that is simple enough so that the par- 
ticle interactions can be described theoretically, but ex- 
hibits the full range of glassy behavior as auto-correlation 
functions with two step relaxation processes showing di- 
vergent relaxation time scales on approaching the glass 
transition. 

In 1984, mode couping theory (MCT), which is a mi- 
croscopic theory for glassy dynamics, was introduced by 
Bengtzelius, Gotze and Sjolander and Leutheusser 0, y] ■ 
It was able to correctly predict many features of the 
complex dynamics of glass-forming liquids and studied 
in the subsequent two decades in great detail by Gotze 
and coworkers; see Refs. [1, Q for detailed reviews. The 
basic version of MCT considers isotropic and homoge- 
neous one-component liquids in three spatial dimensions 
(3D). The only model-dependent input is given by the 
static structure factors of the considered liquid. The 
most prominent prediction of MCT is a dynamic tran- 
sition from a liquid into an ideal nonergodic glassy state 
at some critical temperature or particle density. 

In reality, one-component (monodisperse) systems do 
not serve as good glass-formers since they tend to form 
crystals rather than amorphous solids. Crystallization 
can be suppressed by using polydisperse systems. The 
simplest polydisperse system is a binary mixture. Since 
it is well-known that adding a second component to a one- 
component liquid may strongly influence both its static 
and dynamic properties, Gotze and Voigtmann have 
investigated systematically the glass transition behavior 
of binary hard spheres in 3D. They have found four mix- 
ing effects: (i) small size disparities stabilize the glass, (ii) 
large size disparities stabilize the liquid, increasing the 
concentration of the smaller particles leads to both (iii) 



an increase in the plateau values of the normalized cor- 
relation functions for intermediate times for wave num- 
bers that are not too small and (iv) a slowing down of 
the relaxation of the correlators of the bigger particles 
towards their plateaus. These results qualitatively agree 
with those from dynamic light scatterin g experi ments 0- 
Q and molecular dynamics simulations jlOllTl|. 

Several physical phenomena like equilibrium phase 
transitions strongly depend on the spatial dimensionality 
v. Thus, there naturally appears the question about the 
I?-dependence of the glass transition. From a fundamen- 
tal point of view, there are interesting studies concern- 
ing glass transitions in high dimensions, see for instance 
the recent publications of Schmid and Schilling [T^] or 
Ikeda and Miyazaki [l^ and the references therein. In 
the present paper, we will restrict ourselves to the case 
V = 2. An experimental realization of a model glass- 
former in two dimensions (2D) was presented by Ebert 
et al. [ij]. They consider binary mixtures of super- 
paramagnetic colloidal particles confined at a water-air 
interface which interact via repulsive dipole-potentials. 
The magnetic moments are induced by an external mag- 
netic field perpendicular to the water interface. The re- 
sults for time-dependent correlation functions measured 
by video microscopy clearly exhibit slow glassy dynam- 
ics as found by Konig et al. [l5| and Mazoyer et al. fl6l . 
Computer simulation results of Santen and Krauth [13| 
for polydisperse hard disks in 2D also give evidence for 
the existence of a dynamic glass transition in 2D. Bayer et 
al. [3 have explored the question on the I?-dependence 
of the glass transition by solving the mode-coupling equa- 
tions for a one-component system of hard disks in 2D. 
They have found an ideal glass transition. On a qualita- 
tive level, the results of Bayer et al. for the glass transi- 
tion scenario for monodisperse hard disks in 2D are very 
similar to corresponding MCT results of Franosch et al. 
(l9| for one-component systems of hard spheres in 3D. 

The MCT study of Bayer et al. [Hi was extended to 
binary mixtures of hard disks by Hajnal et al. [l^ and 
also to binary mixtures of dipolar particles in 2D 
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For binary hard disks in 2D the same four mixing ef- 
fects occur as have been reported before by Gotze and 
Voigtmann Q for binary mixtures of hard spheres in 3D. 
Furthermore, it was shown that the glass transition dia- 
gram for binary hard disks in 2D strongly resembles the 
corresponding random close packing diagram. This fact 
is a hint for the applicability of the MCT approximations 
in 2D. 

So far, the MCT results of Hajnal et al. for the 
dynamics of binary hard disks in 2D have not yet been 
tested systematically within the framework of atomistic 
computer simulations. To fill in this gap is the main mo- 
tivation of our contribution. For this purpose we per- 
form Brownian dynamics (BD) simulations for binary 
mixtures of hard disks in 2D. First, we verify the exis- 
tence of the four mixing effects predicted by MCT which 
were briefly described above. Second, we present a quan- 
titative comparison of time-dependent collective density 
correlators from MCT and our BD simulations. 

In order to achieve a self-contained presentation, we 
have organized the paper as follows: in Sects. II and III 
we introduce the correlation functions of central interest 
and we define our model system. In Sect. IV we describe 
our BD simulation techniques. Sect. V contains a brief 
review of the equations and central predictions of MCT. 
We present our results in Sect. VI. We summarize and 
conclude in Sect. VII. 



II. PRELIMINARIES 



particle density of the component a of the liquid where 
ra.iit) is the position of particle i of the component a at 
time t and 5[-] is the Dirac delta distribution. The time- 
dependent density fluctuation of the component a of the 
liquid at the wave vector fc 7^ is given by the spatial 
Fourier transform 



^exp[ifc • rajit)] 



(1) 



of the particle density n°'{r,t). We focus our discussion 
on the matrix of time-dependent partial autocorre- 
lation functions of density fluctuations, which provide a 
statistical description of a multicomponent liquid. For 
i > 0, its components at wave number k are defined by 
the expressions 



= ,1m <i>fW, 



fc > 0, (2) 
(3) 



where 



ITL 



stands for canonical averaging followed by 



carrying out the thermodynamic limit. The zero-time 
value of the correlation matrix defines the normalization 
#(0) = S > 0, the positive-definite static structure fac- 
tor matrix whose elements obey limfe_j.oo S'^^ = XaSa/s- 
Here Sa/s depicts the Kronecker delta and Xa = N^/N 
the particle number concentration of the component a. 



A. Matrix algebra 



III. MODEL SYSTEM 



In the following, wc will make use of the compact math- 
ematical notation introduced in Ref. [2^. Bold symbols 
A, B etc. denote arrays of M x M matrices whose 
components Ak, Bk are labeled by subscript Latin in- 
dices. Their elements A'^^ , B'^^ , also denoted by {A)'^^ , 
(B)^^, are indicated by superscript Greek indices. Ma- 
trix products are defined component-wise, e.g. C = AB 
means Ck = AkBk for all k. A is called positive- (semi- 
)definite, {A ^ 0) A )^ ii this is true for all Ak- For 
discretized model systems where k is restricted to a finite 
number of values, we define the standard scalar product 
{A\B) = EkEo.A'^fyBf where the superscript * 
stands for complex conjugation. The standard norm of 
A is then given by |^| ~ y^{A\A). 



B. Density correlators 



In this work we investigate binary mixtures of hard 
disks in 2D with diameters da which arc distributed 
isotropically and homogeneously with total particle num- 
ber density n. Consisting of "big" {a — b) and "small" 
(a = s) particles with diameters ds < db, the system is 
coupled to a heat bath with thermal energy fc^T and its 
dynamics is governed by Brownian motion. The masses 
TOq and the single-particle short-time diffusion coeffi- 
cients 1?° are set to rus ~ mi, = mo and D'j = DjJ = Dq, 
for simplicity. 

It is well-known, that the thermodynamic equilibrium 
state of the considered model system depends on three 
independent control parameters. Making use of this im- 
plicitness we choose them to be the total 2D packing 
fraction ip = n{Tr/ A){xbdl + Xsd^, the particle number 
concentration Xg =■ Ng/N of the smaller disks, and the 
size ratio S — dg/di,. 



We consider an isotropic and homogeneous classical 
fiuid consisting of M macroscopic components each con- 
taining Na particles of a species a in V spatial di- 
mensions. The total number of particles in the sys- 



tem is then given by = Sa^i-^a- ^'^^ n°'{r,t) 



Si^Ti S['r—ra,i{t)] denote the time-dependent microscopic 



IV. BROWNIAN DYNAMICS SIMULATION 

The basic concept of the algorithm has been described 
in detail in three dimensions in and can easily be 
adapted to the 2D model specified above [l^. We con- 
sider binary mixtures of hard disks with the size ratios of 
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5 ~ dg/df, G {5/7,1/3} with particle number concentra- 
tions Xs e {0.4,0.5,0.6,0.7,0.8}. N = 1000 hard disks 
move in a 2D simulation box of volume V with periodic 
boundary conditions at packing fraction if as defined in 
section lllll After placing the particles on their initial 
positions we provide Gaussian distributed velocities with 
variance (Iwa.iP) = Vq. To propagate the system at time 
t forward in time, we employ a semi-event-driven algo- 
rithm. For every particle, e.g. for particle i of species a 
at the time t, the algorithm determines the possible col- 
lision time Ai"'^ with any other particle. This is easily 
achieved by solving the equation 



where r 



da + dp 



(4) 



ra,i denotes the vector pointing 
from the center of particle i of species a to the cen- 
ter of particle j of species 



de- 



notes the corresponding relative velocity. The smallest 
solution Ai = min{At^^} for all particle pairs deter- 
mines the next event in the algorithm. All particles can 
then be propagated with constant velocity according to 
raA{t + t') = ra,i{t) + Va,xt' for all t' e [0, At]. At time 
t + At, for two colliding particles the elastic binary colli- 
sion laws impose new velocities 



vp,. 



-*api 



fap. 



\ I- 



-*OiB \ -ta.B 



[r- ■ V- )r- ■ . 



(5) 



(6) 



Due to the boundary conditions any particle in the vicin- 
ity of the box-boundary can collide with an image par- 
ticle coming from the other end of the box with the size 

So far. the algorithm described above yields ballistic 
motion. In order to mimic Brownian motion we modify 
it by introducing a thermostat which at every integer- 
multiple of the time TsVo/ds = 0.01 triggers a so-called 
Brownian step. In the Brownian step, all particle veloci- 
ties are freshly drawn from a Gaussian distribution with 
variance rriQVQ/ {ksT) — 2 for all particles. This assures 
that the particles move diffusively with a short-time dif- 
fusion coefficient Do/{vods) = 0.005 on time scales which 
are large compared to tb- 

As the system starts from a cubic lattice it is necessary 
to wait for the system to relax before meaningful station- 
ary averages can be taken. Equilibration was performed 
with Newtonian dynamics (without imposing the Brown- 
ian step) for 10^ time steps in units of dg/vo- We assume 
that the system is equilibrated, when the time-dependent 
correlation functions do not depend on the time origin. 
Correlation functions (with imposing the Brownian time 
step) were measured in a time window of 10^ time steps 
in units of ds/vo which is equivalent to 2551.02 time steps 
in units of dl /Dq . The collective density correlation func- 
tions given by Eq. ^ can directly be calculated from the 
particle positions. 



We selected systems with radius ratios 6 = 5/7 
and S = 1/3. Simulations were performed for Xg G 
{0.1, 0.2, 0.9} and the structure factors and radial 
distribution functions were compared with the Percus 
Yevick results making it possible to exclude the ones with 
crystallization. For (5 = 5/7 and Xs € {0.4, 0.5, 0.6, 0.7} 
and for 6 = 1/3 and x^ G {0.5, 0.6, 0.7, 0.8} we found 
suitable candidates which are still amorphous at high 
packing fractions ip and thus allow us to investigate the 
glassy behavior. 



MODE-COUPLING THEORY 



A. Basic equations 

The mode-coupling theory (MCT) is based on the ex- 
act Zwanzig-Mori equation with a subsequent application 
of the mode-coupling approximations [j| . For Brownian 
dynamics it reads 

r*(t) + S'"i*(i)+ /" dt'mit-t')^{t') =0. (7) 
Jo 

The components of the matrix of microscopic relaxation 
times T shall be approximated by t^^ = 8ap / [k"^ D'^Xa) 
where D° is the single-particle short-time diffusion coef- 
ficient of a tagged particle of species a inside the fluid. 
MGT approximates the memory kernel m{t) by a sym- 
metric bilinear functional 



(8) 



For a multicomponent liquid in 2? > 2 spatial dimensions 
it reads [201 



■D-l 



E 



a' ,P',a" ,p" 



pk+p 

dp / dq 

Ak-p\ 



^yap-,a'P',."P"^a'p'yc."p" (g) 



where the so-called vertices are given by 



Q:^;a (5' ,ol" 13" 



k;p,q 



n Pq aa'a" PP'P" 



CP, _ (fc' +P^~ q^)cf5a, + {k^ -p^+ q^y^.^Scp 



kpq 



(11) 

c^^ denote the direct correlation functions and Ox) = 
27r'^/2/r(2?/2) is the surface of the P-dimensional unit 
sphere. r(a;) is the gamma function, c is related to 5' via 
the Ornstein-Zernike equation (5'"^)^'^ = Sap /xa — nc^^ . 



B. Discretized model 

For practical purposes, we follow Ref. [l^l and dis- 
cretize the wave number fc to a finite, equally spaced grid 
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of K points k = {dd + k)Ak with k = 0,1, . . . , K — 1 and 
< o^; < 1. The integrals in Eq. © are then replaced 
by Riemann sums 

/ dp / dq...^{Akrj2 E •■• (12) 

^I'^-^'l P=o 4-\k-p\ 

and Eq. (jT]) represents a finite number of coupled nonlin- 
ear integro-differential equations. 

The discretized MCT equations shall be applied to the 
binary hard disk mixture in 2D introduced in Sect. IIIII 
For this, we choose K = 250, Ak-di, = 0.3 and 62 = 0.303. 
Static input for the MCT calculations presented in this 
paper is taken from our BD simulations. For the numer- 
ical solution of Eq. ([7]) we use the algorithm described in 
(2^ . Our time grids consist of 256 points and as initial 
time step size we choose 2.5 • 10~^^ units of dl/Do. 

C. Glass transition singularities 

The non-ergodicity parameters (NEPs) are the ele- 
ments of the matrix F = limf^oo^(i)- For the dis- 
cretized model above, it can be proved [2^ that F is 
(with respect to ^) the maximum real, symmetric fixed 
point of the nonlinear map 

I[X] = S - {S-^ + J^[X , X])-\ (13) 

For the binary hard disk model, MCT predicts that the 
liquid-to-glass transition takes place at a locally smooth 
critical surface H which we can represent as the critical 
packing fraction Lp'^{xs, 5) as a function of the concentra- 
tion of the smaller disks and the size ratio. At this surface 
F jumps from (liquid) to some F'^ >- (glass). Quan- 
tities corresponding to critical points shall be indicated 
by a superscript c in the following. General properties of 
% have been discussed in detail in Ref. [2^ . 

Besides F'^, there are further important quantities 
characterizing a generic liquid-to-glass transition point 
which is an singularity according to the classification 
of Arnol'd [2g. Linearization of X*^ around F"^ yields a 
so-called positive linear map [2^ 

C'lY] = 2{S'- - f=):f=[f^ y](S"= - f") (14) 

with ClY] t for aU Y h 0. This map has a non- 
degenerated maximum eigenvalue r = 1 with a corre- 
sponding (right) eigenvector H'^ and a corresponding 
left eigenvector H'^ which is an eigenvector to eigenvalue 
r* = 1 of the adjoint map of C with respect to the scalar 
product defined in Sect. Ill Al These two eigenvectors are 
determined uniquely if we require the normalization 

(H'=|ff'=) {H''\H''{S'' - F=}-iif'=) = 1. (15) 

A further important quantity is the so-called exponent 
parameter 

A'^ = (H'^liS'' - F=}:P"[i^^ H'^]{S'- - F^}) (16) 



whose value determines the exponents in the asymptotic 
scaling-laws (see the next section). These positive expo- 
nents are the critical exponent a obeying the relation 

r2(l-a)/r(l-2a) = A", (17) 

the von Schweidler exponent b satisfying 

r2(i + 6)/r(i + 26) = A^ (18) 

and the exponent 

7 = (a + 6)/(2a6) (19) 

describing the divergence of the time scale for the final 
relaxation of ^{t) to upon increasing ip towards its 
critical value ip"^. 

D. Asymptotic scaling-laws 

Close to the liquid-glass transition, MCT makes uni- 
versal predictions for the relaxation behavior of $(i) 
which can be studied in the framework of asymptotic 
expansions. For the following, let us fix Xg and 6 to some 
specific value and define the distance parameter 

£ = (^-^^)/^^ (20) 

Finally, we introduce the separation parameter a{s) 
which is a linear function of e. It follows from 

aiip) = (Hulls'- - F'}S^~^{ST[F^,F^]{S - F^) 
-5'":P'=[F'=,F'=](S"=-F=)}) (21) 

by expanding around 95^ up to linear order in e 0, [23| • 

1. The first scaling-law regime 

For small separation parameters, #(t) develops power- 
law dynamics located around F'^ which is also called the 
/3-relaxation process. For times within the so-called first 
scaling-law regime defined by \^{t) — F'^\ <^ 1, also called 
the /3-scaling regime, there holds the factorization theo- 
rem 

^{t)^ F^ = H^git)+Oi\a\) (22) 

with the /3-correlator Q{t) ~ 0{\a\^l'^) obeying the equa- 
tion of motion 

o + \^g''{t) = ^ ^* dt'g(t - t')g{t') (23) 

with the divergent initial condition 

g(t ^ 0) = (tmY\ (24) 

The e-independent time scale io has to be matched to 
the full solution of Eq. ([7]) at e = since in this case the 
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power law occurring in Eq. (|24l) is a special solution of 
Eq. which describes the relaxation of ^{t) towards 
F'^-. It is easy to verify that G(t) obeys the scaling-law 

g{t, a ^ 0) - W\^^^-g{i = t/t,,a = ±1), (25) 

U=to\a\-i^. (26) 

The master function Q{t) obeys Eq. with the re- 
placement a I— )■ cr/|cr| = ±1 and the initial condition 
g{i ^ 0) = (i)~°. Eq. ^ can be solved by asymptotic 
series expansions j28| . One finds the leading long-time 
asymptotes 

g{t ^oo, a >0)^^C7/{1-X-), (27) 

g{t^ oo, a <0) = -{t/T)\ (28) 

Eq. (P7)) describes the asymptotic behavior of F in the 
glassy regime close to (p'^. Eq. (|28p is referred to as the 
von Schweidlcr law and describes the initial part of the 
relaxation of $(t) from to 0. Eqs. and (^5]) imply 

T-ftolfrr^ (29) 

with the £r-independent constant f which demonstrates 
that the arrest of ^{t) to F'^ at ip = ip'^ is caused by a 
power-law divergence of the times scale for the inset of 
the relaxation of $(i) from F'^ to upon increasing p 
towards its critical value . 

2. The second scaling-law regime 

Now we turn to the relaxation process of $(<) from 
to within the liquid regime which is also called the 
a-relaxation process. Considering the limits e — >■ 0~ and 
i — oo with fixed t = t/r, one arrives at the so-called 
a-scaling law [1, [l^ 

*(t) = *^(t)+0(H) (30) 

where the e-independent master function #'^(f) obeys the 
equation of motion 

d /■* 

(5'c)-i$'=(f) ^ rh=(f)5^ - / d?m'=(f- ?)*"(?) 

di Jo 

(31) 

with the memory kernel 

m%i) = J^'^l^'^ii), *=(?)] (32) 
and the von Schweidler law 

*^(t ^ 0) = F'^ - {t)''H'' (33) 

as initial condition. Eq. ([50)1 is also called superposi- 
tion principle due to the following implication: for times 
within the so-called second scaling-law regime given by 
i ^ T, density correlators $^'^(i) corresponding to dif- 
ferent values of e collapse onto master curves (#'^)^^(i) 
when they are plotted as functions of t = t/r. 



E. Mixing effects 

As already briefly mentioned in the introduction, MCT 
predicts for binary hard disks in 2D the existence of four 
mixing effects [20]. Let us recapitulate these effects in 
more detail. 

(i) For small size disparities the glassy regime is en- 
hanced. For 0.65 < (5 < 1 it is (p'^{xs,S) < if 
< < 1, where ip^ denotes the critical pack- 
ing fraction for monodisperse hard disks in 2D. For 
S = 5/7, for instance, ip^ixg) develops a single min- 
imum located at Xs ~ 0.5. 

(ii) For larger size disparities the liquid state is sta- 
bilized, i.e. for S < 0.37 it is p''{xs,S) > (pQ if 
< Xs < I- This effect is also called plasticization. 
For (5 = 1/3, for instance, (p'^{xs) develops a single 
maximum located at Xs ~ 0.85. 

(iii) Upon increasing the concentration Xs of the smaller 
particles the NEPs, and thus also the plateau values 
of the normalized correlation functions ^'i^"{t)/S^" 
in the liquid regime for intermediate times, increase 
for not too small k and all Xs- 

(iv) Starting with = and increasing the concen- 
tration Xs of the smaller particles leads for not 
too large fc to a slowing down of the relaxation 
of the normalized correlators '^'^ {t) / S'j}' of the big 
particles towards their plateaus in the sense that 
the ^f'{t)/Sl'' versus \ogiQ{t) curve becomes flat- 
ter upon increasing Xg. 

The amplitude of the predicted effects (i) and (ii) is quite 
small, the total variation of (p"^ {xg , S) is of the order of 
1%. Nevertheless, these small changes in ip'^ may have 
a strong influence on time-dependent quantities which 
are accessible to our BD simulations. Let us fix some 
S and ip such that for all < x^ < 1 the condition 
< [p}''{xs,S) - </j] < 1 is satisfied. Eqs. ^ and §UI 
predict then a strong variation in the time scales r for 
the a-relaxation of the correlators ^'^^{t). It has been 
demonstrated [20j that an occurring minimum (maxi- 
mum) in ip''(xs) is then directly reflected by an occurring 
maximum (minimum) in r. Using this information, we 
demonstrate in the following that the mixing effects (i)- 
(iv) predicted by MCT are indeed observable in our BD 
simulation data. 



VI. RESULTS AND DISCUSSION 

We choose to use a Brownian dynamics simulation for 
our comparison with MCT results for the following rea- 
son: the behavior of the relaxation times of the system on 
approaching the glass transition point, the details of the 
a-relaxation, as well as the NEPs are independent of the 
microscopic dynamics of supercooled liquids, as Gleim et. 
al. showed in [2^. Nevertheless, the relaxation onto the 
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Kohlrausch laws, 



(34) 



FIG. 1: (Color online) Normalized critical non-ergodicity pa- 
rameters of the simulated collective density correlators ^'k'{t) 
of the big particles, extracted from Kohlrausch fits for the size 
ratio 5 — 5/7 and the packing fraction ip — 0.79. The num- 
ber concentrations of small particles Xs vary as labeled in the 
legend. The solid red line gives the MCT results calculated 
with a simulated structure factor input at flicT ~ 0.6920 
and Xs — 0.5. 



0.4 
□ xl = 0.5 
o ^, = 0.6 
A x^ = OJ 

0.5 (MCT) 




FIG. 2: (Color online) Normalized critical non-ergodicity pa- 
rameters of the simulated collective density correlators (t) 
of the small particles, extracted from Kohlrausch fits for the 
size ratio S — 5/7 and the packing fraction ip = 0.79. The 
number concentrations of small particles Xa are the same as 
in Fig. [T] The solid red line shows the MCT results calculated 
with simulated structure factors as in Fig. [1] 



plateaus is expected to be described better within MCT 
for Brownian dynamics (29| . 



A. Statics 

The first point we want to address are the non- 
ergodicity parameters. A common description of 
the a-relaxation is in terms of stretched exponential 



with the stretching exponent , a relaxation time scale 
f^" and the amplitude A'^". For structural relaxation 
in equilibrium system /3^" < 1 is required. The a- 
master function from MCT equation (|3ip is different 
from the Kohlrausch form, however the theory predicts 
that for large wave numbers the two functional forms 
become identical and /3^" — b [sO]. The Kohlrausch 
amplitude ^4^" provides an estimate for the MCT NEPs 
(i^"^)^". Since the a-process starts below this plateau 
value, < (i^"^)^" should hold. However, in practice 
the separation of the a-process from the /^-relaxation is 
not clear enough to fulfill this prediction. 

Kohlrausch fits are hindered by some subtle problems: 
lacking a clear separation of the a-process, the fit pa- 
rameters inclose a dependence on the fit range. A priori 
it is unclear how to choose the optimal fit range, as for 
very long times one expects the relaxation to become 
(non stretched) exponential again, and for short times, 
deviations stemming from the /3-relaxation hamper the 
choice. The fit range was fixed, so that the parameters 
only exhibit the weakest (the region were they are almost 
constant) dependence on the boundaries. This procedure 
leads to t e [7.65,2551.02] for S = 5/7 and if = 0.79 
with the various Xs- For 6 = 1/3 and (p ~ 0.81 we used 
t e [1.39,555.56] for the various Xg- 

Fig. [1] shows approximate values for the normalized 
critical NEPs {F^)f/{S'=)f for the big disks at S = 5/7 
and different values for Xg , extracted from our BD simu- 
lation data via Kohlrausch-fits. Corresponding results for 
the smaller disks are shown in Fig. [51 In both Fig. [1] and 
Fig. [21 we have also included critical NEPs for Xs = 0.5, 
calculated via MCT with BD-simulated structure factors 
as input, with MCT yielding a critical packing fraction of 
ftiCT ~ 0.6920 for the simulated input. On a qualitative 
level, our MCT results are in good agreement with our 
BD simulation results. For the big particles the relation 
j^aa ^jpc-^aa jg fulfilled for all kdb except for some 
outliers. The same holds for the small particles but for 
kdi, < 5 the Kohlrausch fit yields smaller estimations for 
the NEPs. In both Fig. [Hand Fig. [21 we observe a slight 
increase in the BD simulation results for the NEPs upon 
increasing Xs which is on a qualitative level in agreement 
with previous MCT results @,[23]. 

Fig. [31 shows approximate values for the normalized 
critical NEPs {F'')^ / {S^f^ for the big disks at <5 = 1/3 
and different values for Xs , extracted from our BD simu- 
lation data via Kohlrausch-fits. Corresponding results for 
the NEPs of the smaller disks are shown in Fig.[31 In both 
Fig.[31and Fig.[31 we have also included critical NEPs for 
Xs = 0.5, calculated via MCT with BD-simulated struc- 
ture factors as input, giving a critical packing fraction of 
'fi'iiCT ~ 0-6991 for the simulated input. For the present 
value of 6, we observe that our MCT calculations yield 
systematically smaller values for the NEPs, compared to 
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FIG. 3: (Color online) Normalized critical non-ergodicity pa- 
rameters of the simulated collective density correlators ^'k'{t) 
of the big particles, extracted from Kohlrausch fits for the size 
ratio S = 1/3 and the packing fraction ip — 0.81. The vari- 
ous number concentrations Xs are as depicted in the legend. 
MCT results using simulated structure factors as input with 
vliCT = 0.6991 are shown as black solid line for Xs = 0.5. 




FIG. 4: (Color online) Normalized critical non-ergodicity pa- 
rameters of the simulated collective density correlators {t) 
of the small particles, extracted from Kohlrausch fits for the 
size ratio 5 = 1/3 and the packing fraction = 0.81. The 
color and symbol coding for the different Xa is the same as in 
Fig. El The solid black line shows the MCT results calculated 
with simulated structure factors as used in Fig. [31 



our BD simulations results. The underestimation of the 
NEPs may be attributed to the underestimation of tp'^: 
MCT predicts arrest at lower densities, but the NEPs 
may increase with density as the denser glass is stifFer 
with respect to density fluctuations. The BD simulation 
results in both Fig. [3] and Fig. S] indicate a systematic 
increase in the NEPs upon increasing Xs which is, as 
expected, more strongly pronounced than for the case 
6 = 5/7 which is on a qualitative level in agreement with 
previous MCT results flU^. 




FIG. 5: (Color online) Critical amplitudes obtained by 
Eq. H35|l for the big and small particles with kodb = 2.5. 
The data were extracted from the collective correlators at 
(f = 0.79, 5 = 5/7 and Xs = 0.5. Solid black and dashed red 
lines depict the MCT results obtained with the same simu- 
lated structure factors as in Fig. [H and Fig. (21 
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FIG. 6: (Color online) Critical amplitudes obtained by 
Eq. (|35p for 5=1/3 and Xs = 0.5 with kodb = 5.37 from 
the simulated correlators a,t ip = 0.81, 5 = 1/3 and Xs = 0.5. 
Solid black and dashed red lines depict the MCT results ob- 
tained with the same simulated structure factors as used in 
Fig. [Hand Fig.H 



We can conclude here, that at least on a qualitative 
level, our BD simulation results confirm the Xs and i5 
dependences of the NEPs predicted by MCT In 
particular, we have clearly verified the existence of mixing 
effect (iii). 

Let us investigate the so-called critical amplitude as 
a further interesting static quantity. Fig. [5] shows ap- 
proximate values for the normalized critical amplitudes 
(ff'=)^"/(S"=)^" for both the big and the small disks at 
6 = 5/7 and Xg = 0.5, extracted from our BD simu- 
lations. Corresponding results for 6 = 1/3 are shown in 
Fig. [5| In order to determine the critical amplitudes from 
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FIG. 7: (Color online) Simulated normalized collective corre- 
lation functions of the big particles for the size ratio 5 = 5/7 at 
if = 0.79 and kdb = 8.5 for varying Xs G {0.4, 0.5, 0.6, 0.7} as 
labeled in the legend. Solid lines show examples of Kohlrausch 
fits to the Xs = 0.4 correlator (black) and the Xa = 0.7 corre- 
lator (blue). 



FIG. 8: (Color online) Simulated normalized collective corre- 
lation functions of the small particles for the size ratio 5 = 5/7 
at = 0.79 and kdt = 8.5 for varying a::^. The color and sym- 
bol coding is the same as in Fig. [T] Solid lines show examples 
of Kohlrausch fits to the Xs = 0.4 correlator (black) and the 
Xs = 0.7 correlator (blue). 



the BD simulation data one can define the function [3l| 

(35) 

with tj chosen in the /3-scaling regime vifhich in our case 
is determined to tj e [0.47449,8.2882] for 6 = 5/7 and 
tj € [0.10334, 1.8051] for S = 1/3. The last equality fol- 
lows from Eq. (|22|) and thus allows us to extract the criti- 
cal amplitudes {H'^)'^^ up to a factor (iJ^)^^^. Beside the 
numerical uncertainty at low fcdt,, the simulation data for 
i5 = 1/3 are in a good agreement with our corresponding 
results from MCT using BD-simulated structure factors 
as input. For the case ^ = 5/7 we observe larger discrep- 
ancies for all fcdf). 



B. Dynamics 

In this section the most important quantities for 
MCT, the collective density correlators will be discussed. 
Fig. [7] shows normalized collective density correlators 
^^k{t)/S^ from our BD simulations for the big parti- 
cles of binary hard disk mixtures in 2D &i 5 ~ 5/7, 
If = 0.79 and kdi, = 8.5 for different concentrations Xs 
of the smaller disks. Similar results for ^^k{t)/S^'' are 
shown in Fig. [S] Focusing on the data for ^^]^{t)/S^ in 
Fig. in these data exhibit the same three mixing effects as 
the ones from MCT shown in Fig. 6 in Ref. [20], namely 
(iii) an increase in the plateau values accompanied by (iv) 
a slowing down of the relaxation towards these plateaus 
and (i) an additional slowing down of the a-relaxation 
process upon increasing x^- However, there are some de- 
viations: from Fig. 5 in Ref. [13] we would expect that the 



slowest a-relaxation process occurs at Xg = 0.5. The sim- 
ulation data, however, exhibit the slowest a-relaxation at 
the highest investigated value Xg = 0.7. Unfortunately 
the simulation systems at higher Xs are subject to crys- 
tallization which makes them unsuitable for the MCT 
comparison, rendering this region 'unaccessible' for the 
simulations. A possible source for the discrepancy could 
be, that in Ref. [13] an approximate theory, the Percus- 
Yevick structure factor was used as input to the MCT 
calculations. Furthermore, for 5 = 5/7, MCT predicts 
0.686 0.6920 for all x^. The simulation data, 

however, imply ^%i„^ ^ 0.79 which means that MCT un- 
derestimates the critical packing fraction by about 15%. 
The Pcrcus-Yevick approximation contributes to this un- 
derestimation but even with simulated structure factors 
as input, MCT underestimates the critical packing frac- 
tion. In Fig. [H] we basically find the same scenario as in 
Fig. [3 although a bit less pronounced, as expected from 

MCT Hi]. 

Fig. [9] shows the normalized collective density corre- 
lators <^f{t)/Sl'' from our BD simulations at (5 = 1/3, 
if = 0.81 and kd^ = 9.0 for different concentrations Xs 
of the smaller disks. Similar results for ^'ff{t)/S^'' are 
shown in Fig. (TO] First we have a closer look on the 
data for $f (t)/5'f in Fig.lH On a qualitative level, the 
data for the three lowest values for Xs are fully consis- 
tent with all MCT results in both Fig. 5 in Ref. ^ 
and Fig. 7 in Ref. [20j . Upon increasing Xg, the simu- 
lation data exhibit the mixing effect (iii) an increase in 
the plateau values accompanied by (iv) a slowing down 
of the relaxation towards these plateaus. In addition to 
that, increasing Xg from 0.5 to 0.7 leads to (ii) a speeding 
up of the a-relaxation. Thus, the three correlators cor- 
responding to the lowest values for Xg exhibit a pair-wise 
crossing. A further increase in Xg to 0.8 leads again to a 
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FIG. 9: (Color online) Simulated normalized collective corre- 
lation functions of the big particles for the size ratio (5 = 1/3 at 
V9 = 0.81 and kdt, = 9.0 for varying a;^ G {0.5, 0.6, 0.7, 0.8} as 
labeled in the legend. Solid lines show exemplary Kohlrausch 
fits to the Xa = 0.5 correlator (black) and the Xs = 0.8 corre- 
lator (blue). 




FIG. 10: (Color online) Simulated normalized collective corre- 
lation functions of the small particles for the size ratio 5 = 1/3 
at = 0.81 and kdi, = 9.0 for varying Xs- The color and sym- 
bol coding is the same as in Fig. Solid lines show exem- 
plary Kohlrausch fits to the a::^ ~ 0.5 correlator (black) and 
the Xa = 0.8 correlator (blue). 



slowing down of the a-relaxation process, although from 
the MCT results in Fig. 5 in Ref. [l^l we would expect 
the fastest a-rclaxation at Xg = 0.85. For the small parti- 
cles we observe the similar albeit less pronounced effects, 
see Fig. [lOl 

The fact that the simulation data at (5 = 5/7 and f = 
0.79 and the ones a.t S — 1/3 and (p = 0.81 show very sim- 
ilar a-relaxation times is on a qualitative level consistent 
with the MCT result ip''{xs,S = 5/7) < f''{xs,S = 1/3) 
for < < 1. Let us conclude here with the statement 
that, at least on a qualitative level, the four mixing ef- 
fects predicted by MCT for the binary hard disk model in 



2D [20[ are also observable in our computer simulations 
which supports the quality of MCT in 2D. 



VII. THE GLASS TRANSITION OF A 
SELECTED MIXTURE 

In this section we select one of the systems {S = 5/7, 
Xg = 0.5) from the preceeding sections and perform 
a more profound MCT analysis including asymptotic 
checks and determining the MCT glass transition point. 
This will complete the MCT analysis of the system al- 
ready discussed under shear in Ref. [2^ with respect to 
the quiescent state. 

In analogy to Eq. ([55)1 it is possible to test another pre- 
diction of MCT. In order to investigate the factorization 
theorem given by Eq. (1221) . we consider the function (ssj 



$f (i')-1>f (i") 



(36) 



with fixed times t' < t" to be chosen appropriately from 
the /3-scaling regime. Then Eq. predicts 



G{t)~G{t') 
G[t') - Git") 



-0(\a\ 



(37) 



not to be dependent on wave number and particle index, 
to leading order in the separation parameter a. Thus it 
must be possible to fix two times t' and t" uniquely so 
that superimposing X^^{t) for different k yields a win- 
dow in which all X^^ (t) collapse. An advantage of this 
procedure is that the critical amplitude drops out and 
doesn't need to be fitted. 

Fig. [TT] shows our BD simulation results for X^"(t) 
at ^ = 0.79, 5 = 5/7, x^ = 0.5, t'Do/dl = 0.7648 
and t"Do/dl = 9.117 for different wave numbers. In- 
deed, within the numerical accuracy of our simulations, 
for both a — b and a = s the data for different kdh col- 
lapse onto each other within a time window of about two 
decades, similar to previous results for binary Lennard- 
Jones mixtures in 3D [s^ . A more sensitive test of MCT 
asymptotics is the so-called ordering rule. As in the next- 
to-leading order corrections to the factorization theorem 
the same fc-dependent correction amplitudes appear, the 
deviations before the collapse regime must be in the same 
direction as after the collapse window. Hence correlators 
entering the collapse region in a certain order when num- 
bered from top to bottom should leave the collapse win- 
dow in exactly that ordering [l^. Fig. [11] is confirming 
that prediction. 

Now, we test the validity of the a-scaling law given 
by Eq. ((30)) . According to that equation, plotting the 
correlators as a function oi t/r makes the data collapse 
for long times on a master curve on approaching ^ Lp'^ 
from the liquid. To determine a dimensionless relaxation 
time r/r* oc \i7\~^ from the simulation, the correlators 
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FIG. 11: (Color online) Functions X^°'{t) calculated from 
Eq. (|36p with the simulated correlators at 93 = 0.79, S = 5/7 
and Xs = 0.5 by fixing t'Do/dl = 0.7648 and t"Do/dl = 9.117 
for the big and the small particles. 
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FIG. 12: (Color online) Rescaled correlators for tp < 0.79, 
S = 5/7 and Xs ~ 0.5 to collapse on one a-master function at 
long times. The rescale times are independent on kdb 



corresponding to Lp < 0.79 have been shifted along the 
log2Q(i)-axis to coincide in the final decay with the one 
at the highest packing fraction ip = 0.79. The highest 
packing fraction defines Tip^o.Tg/'^* = 1 a-nd yields our 
best approximation for the a-master function. Checking 
that t/t* is independent on kdf, we can validate the a- 
scaling. We have chosen the structure factor peak for the 
determination of t/t* , as the strength of the a-process is 
maximal here and thus a separation from the /3-process 
can be achieved. Fig. [12] shows the exemplary result of 
the a-scaling for two different wave numbers. The data 
clearly exhibit the two-step relaxation pattern of glass- 
forming liquids with increasingly stretched plateaus upon 
increasing the packing fraction. The shifted correlators 
approach an a-master curve with the highest densities 
collapsing over almost three decades in time. 

After having checked a few asymptotic results we now 




FIG. 13: (Color online) Normalized total collective correlators 
from BD simulations for S = 5/7, Xa = 0.5 and kdb = 8.5. The 
solid lines show corresponding MCT results where the packing 
fractions have been mapped according to Eq. (|38p . The inset 
shows the a-relaxation timescales and the fitted power law 
t/t* = A\ip - ipsiml^'' with 7 = 2.4969 (fixed from MCT), 
A = 1.633 • 10"'^ and (fi%^ = 0.7948. See text for details. 
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FIG. 14: (Color online) Normalized collective correlators from 
BD simulations for S — 5/7, Xs = 0.5 and kdb = 8.5. The 
solid lines show corresponding MCT results where the packing 
fractions have been mapped according to Eq. (|38[) . 



present full numerical MCT calculations. As a prelimi- 
nary we use BD-simulated static structure factors as in- 
put for MCT to calculate the critical packing fraction 
^pIjct - 0.6920 and the exponent 7 = 2.4969. The ob- 
vious mismatch in ip'^ necessitates a comparison, at cor- 
responding separation from the transition point. That 
entails matching the separation parameter a which is not 
easy to obtain from the simulation, but a peculiarity of 
our system helps us to circumvent the problem. In Fig.fT^ 
we show BD simulation results for the normalized total 
collective correlators X]a/3 '^fe''(^)/ '^fc'' at (5 = 5/7, 
Xs = 0.5 and kdb = 8-5 for different packing fractions 
(p within the liquid regime close to vitrification. The in- 
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FIG. 15: (Color online) Normalized collective correlators from 
BD simulations for S = 5/7, Xs = 0.5 and kdb = 6.4. The 
solid lines show corresponding MCT results where the packing 
fractions have been mapped according to Eq. (|38|l . 




FIG. 16: (Color online) Normalized collective correlators from 
BD simulations for S = 5/7, Xs = 0.5 and kdb ~ 13.0. The 
solid lines show corresponding MCT results where the packing 
fractions have been mapped according to Eq. (|38|l . 

set in Fig. [131 shows our results for t/t* obtained from 
the shifting process seen in Fig. [121 The straight Hne 
shows the result from a power- law fit t If ^ Vsim\~'^ 
with fixed 7 = 2.4969 which yields the extrapolated value 
fsirn — 0.79481 for the critical packing fraction which is 
approximately 15% larger than the value predicted by 
MCT. For - ip%^\ < 0.01 the <y9-dependence of t/t* 
is excellently described by the MCT exponent 7. We 
can conclude: although MCT underestimates the critical 
packing fraction, it nevertheless describes very well the 
(/^-dependence of the a-relaxation process in the liquid 
regime close to the glass transition, which again supports 
the quahty of MCT in 2D. 

Using this information, we are able to present a quan- 
titative comparison of time-dependent correlation func- 
tions from MCT to those from our BD simulations. For 



this purpose we have to take into account that MCT 
overestimates glass formation. It is well known that this 
results in predicting the glass transition at a too low crit- 
ical packing fraction ip'^. Hence the relevant parameter 
when comparing MCT and simulation results is the sep- 
aration parameter a which depends linearly on the dis- 
tance parameter via a = Ce (sec Section fypp . As the 
constant C is evaluated at the critical packing fraction 
with the corresponding structurefactor it is reasonable to 
assume that MCT doesn't yield the same prefactor as the 
simulation. For instance Flenner and Szamel found that 
both prefactors differ in a 3D binary Lennard Jones mix- 
ture [3JI . In order to construct a mapping of the packing 
fractions fsim used in our BD simulation onto some ap- 
propriate ones fMCT to be used for the corresponding 
MCT calculations we postulate that the separation pa- 
rameters for both systems must be equal. This leads us 
to the Ansatz 

£MCT = {Csim/CMCT) (38) 

with some appropriately chosen constant A, which in our 
special case can be found empirically as A = 1. With 
this, all input parameters for our MCT equations are 
uniquely determined. The solid lines in Fig. [131 represent 
our MCT results for the normalized total collective cor- 
relation functions corresponding to the shown BD sim- 
ulation data in Fig. [T21 Corresponding results for the 
normalized partial correlators {t) j S^" arc shown in 
Fig. [TH We observe that for the chosen wave number 
MCT tends to underestimate the correlation functions in 
the transient time regime tDo/dl = 2.5- 10^"^ and overes- 
timates the plateau values at times within the /3-scaling 
regime. Beside these quantitative deviations, MCT de- 
scribes very well the qualitative t and ip dependences 
of the BD-simulated correlation functions for (p > 0.78. 
Note especially that in this parameter regime MCT de- 
scribes the final part of the a-relaxation process also on 
the quantitative level correctly. 

Of course, the grade of quantitative compliance of our 
BD simulation and MCT results is also dependent on 
the wave number. We demonstrate this in Fig. \W\ for 
^rW/^k" at kdb = 6.4 and in Fig. [HI for $rW/'5fc" 
at kdb = 13.0. For instance, MCT strongly overesti- 
mates the plateau values and also the a-relaxation times 
for ^t{t)/Sf at both kdb = 6.4 and kdb = 13.0. For 
^rlOZ-^T at kdb = 6.4, on the other hand, MCT 
underestimates both the plateaus and the a-relaxation 
times. For $^(07-5^ at kdb = 13.0 MCT slightly over- 
estimates the plateaus, but excellently describes the a- 
relaxation processes. All these data indicate a connection 
between the (over)under-estimation of the plateaus and 
the (over)under-estimation of the a-relaxation times by 
MCT. 

For completeness we show the Kohlrausch stretching 
exponent in Fig. [ITl for the small and big particles. For 
large kdb the exponent converges as /3J?" — )• 0.6 ± 0.05. 
Unfortunately noisy data prevents a more precise deter- 
mination of the high-fcc?b limit. Nevertheless the high-kdb 
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FIG. 17: (Color online) The Kohlrausch stretching exponent 
/3fc" for the big (black circles) and the small (red squares) 
particles at 5 = 5/7, Xs = 0.5 and (p = 0.79. The MOT 
von Schweidler exponent b — 0.5571 is indicated by the black 
dashed line. 
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FIG. 18: (Color online) Three correlators above the MCT 
glass transition point at 5 = 5/7, Xs = 0.5, ip = 0.81, (fi = 0.80 
and — 0.795 together with the highest liquid density at 
= 0.79. The correlators violate the scaling property shown 
in Fig. 1121 The plateau values increase on going deeper into 
the glass 



limit is in good accordance with MCT's von Schweidler 
exponent, Eq. calculated to b = 0.571 with simu- 

lated structure factors. 



VIII. BEYOND THE MCT GLASS TRANSITION 

Fig. [THl shows an example of three correlators with 
packing fractions above the MCT glass transition point. 
Clearly one can see that the plateaus of the correlators 
are increasing with tp. This leads to intersecting cor- 
relators and proves that the correlators do not collapse 
on a master function anymore when they are rescaled 



above ip 



> 



0.7948 show relaxation processes, not 



captured in MCT which make the correlators beyond the 
MCT glass transition point decay to zero for long times. 

The correlator at ip = 0.81 does not show a relax- 
ation within the time window accessible in our BD sim- 
ulations. In the time window accessible by Newtonian 
dynamics a relaxation time could be extrapolated by fit- 
ting a Kohlrausch function to the coherent correlator at 
kdb = 8.4. These correlators were obtained by equilibrat- 
ing 150 independent systems with Newtonian dynamics 
up to 20% of the structural relaxation time measured 
with the Kohlrausch fit. After having performed the 
Newtonian equilibration we do not see any dependence 
on the time origin for correlation functions in the time 
window accessible by Brownian dynamics. Especially the 
plateau values do not change any more. However the on- 
set of a final decay can be assumed, as the curves show 
a slight curvature for large times. 

At this point we want to stress that even though the 
decay to zero is not included in MCT, MCT still gives 
an explanation for the rise of the NEPs above the glass 
transition, see also Ref. 



IX. SUMMARY AND CONCLUSIONS 

We have performed BD simulations for binary mixtures 
of hard disks in 2D in order to test systematically the 
predictions of MCT for this specific model system. Such 
a systematic test comparing collective density correlation 
functions from both approaches has not been carried out 
before for a 2D model system. 

Our main result is that MCT seems to be capable to 
capture many qualitative features of the glass transition 
behavior of binary hard disk mixtures in 2D. Particu- 
larly, all four mixing effects predicted by MCT for bi- 
nary hard disks in 2D are indeed observable in our BD 
simulations. Furthermore, we have demonstrated the va- 
lidity of the factorization theorem given by Eq. and 
the a-scaling law given by Eqs. (|29)) and (pO)) . We have 
found that the MCT exponent 7 = 2.4969 describes ex- 
cellently the dependence of the a-relaxation process on 
the packing fraction ip in the liquid regime close to the 
glass transition. Furthermore, MCT describes very well 
the qualitative t and ip dependences of the BD-simulated 
collective correlation functions close to the glass tran- 
sition if one rescales the packing fraction according to 
Eq. ([55]) . All these facts strongly support the quality of 
MCT in 2D. The range of validity of the asymptotic scal- 
ing law is very similar to previous results for hard spheres 
in 3D. 

Going beyond the extrapolated glass transition at 
'Psim ^ 0.7948 our BD simulation correlators loose the 
a-scaling property which results in intersecting correla- 
tors when they are rescaled by their relaxation times. 

On the quantitative level, we have found some discrep- 
ancies between MCT and our BD simulation results. For 



13 



instance, MCT underestimates the absolute value for the 
critical packing fraction for vitrification by about 15%, 
similar to previous results for hard spheres in 3D. It seems 
to be possible to suppress the influence of this error on 
time-dependent quantities by introducing a linear map 
according to Eq. psp , at least for packing fractions close 
to the glass transition. We have found quantitative devi- 
ations of the BD simulation results for the NEPs and the 
critical amplitudes from the corresponding MCT results. 
For a better understanding of these deviations it would be 
necessary to investigate the leading order corrections to 
the asymptotic scaling laws presented in Sect. IV D] This, 
however, would go far beyond the scope of our present 
study. 

In the present paper we have verified mixing scenarios 
predicted by MCT for binary mixtures of hard disks in 



2D. Recently, mixing effects on the glass transition have 
also been investigated systematically in the framework 
of MCT for binary mixtures of dipolar point particles in 
2D [2l[. As a project for the future it would be worth 
to also test the obtained MCT results systematically by 
BD simulations. 
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